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1 Introduction 



The physics of the mixing and decays of B-mesons is essential for a determination of unknown 
CKM-matrix elements and thus for our understanding of the violation of CP-symmetry in 
Nature. It is also still promising for the discovery of physics beyond the standard model of 
particle physics. Unfortunately, many of the experimental observations can only be related 
to the standard model parameters if transition matrix elements of the effective weak Hamil- 
tonian are known. These matrix elements between hadron states are only computable in a 
fully non-perturbative framework. They provide a strong motivation to study B-physics in 
lattice QCD. However, as the mass of the b-quark is larger than the affordable inverse lattice 
spacing in Monte Carlo simulations of lattice QCD, this quark escapes a direct treatment as 
a relativistic particle. Therefore, effective theories for the b-quark are being developed and 
used to compute the matrix elements in question [1,2]. 

The first — and very promising — effective theory that was suggested is the Heavy 
Quark Effective Theory (HQET) [3, 4]. Like others, it is afflicted by a problem which 
remained unsolved so far: in general its parameters (the coefficients of the terms in the 
Lagrangian) themselves have to be determined non-perturbatively, as briefly explained in 
Sect. l2~2l In other words, the theory has to be renormalized non-perturbatively [5]. This 
fact is simply due to the mixing of operators of different dimensions in the Lagrangian, re- 
quiring fine-tuning of their coefficients. If they were determined only perturbatively (in the 
QCD coupling), the continuum limit of the theory would not exist. 

The issue is already present in the determination of the b-quark mass in the static 
approximation, i.e. in the lowest order of the effective theory. In [6] a strategy was introduced 
and successfully applied to this problem for the first time, and a general framework for a 
non-perturbative renormalization of HQET was sketched in [7]. The basic idea, illustrated 
in Fig. ^ is easily explained. 

In a finite volume of linear extent Lq = 0(0.2 fm), one may realize lattices with am^ <C 1 
such that the b-quark can be treated as a standard relativistic fermion. At the same time 
the energy scale 1/Lq = 0(1 GeV) is still significantly below and HQET applies quanti- 
tatively. Computing the same suitable observables in both theories relates the parameters of 
HQET to those of QCD. Then one moves, by an iterative procedure that we can still leave 
unspecified here, to larger and larger volumes and computes HQET observables. This yields 
the connection to a physically large volume (of linear extent 0(2 fm)), where eventually the 
desired matrix elements are accessible. 

Since in this way the parameters of HQET are determined from those of QCD, the 
predictive power of QCD is transfered to HQET. In addition to solving the renormalization 
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Figure 1: Relating experimental observables to properly renormalized HQET. r sta t is a renormalized 
quantity in HQET and cr m (p 2 (L)) connects r sta t(£) and r sta t(2L); their exact definitions will be given in the 
course of this paper. In the chosen example, the experimental observable is the mass of the B-meson. 
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problem of the effective theory, one also eliminates the usual need to determine more and 
more parameters of the theory from experiment as the effective theory is considered to higher 
and higher order. 

Although related, the strategy we propose here is not to be confused with the one for 
the computation of the running of the coupling and renormalization group invariant matrix 
elements as first suggested by Liischer, Weisz and Wolff [8] and then developed by the 
ALPHA Collaboration. We will discuss the difference in Sect. 15.2.41 

In this paper we define the effective theory in detail, discussing in particular its renor- 
malization properties (Sect. [2j). We then explain the matching between QCD and HQET 
(Sect. |3J) as well as the finite-size strategy (Sect.|IJ in the general case. Sect. El provides two 
examples of applications using the effective theory at the lowest order in the inverse b-quark 
mass. The first one is the computation of the quark mass, where numerical results illustrate 
that indeed the power-law divergence can be subtracted non-perturbatively, retaining a very 
good precision for the final physical number. The second one, devoted to the B-meson decay 
constant, has not yet been applied numerically but is a useful and simple example to help 
in understanding our method. In Sect. El we discuss the potential of our approach as well as 
the expected uncertainty due to the use of a finite order in the HQET expansion. 

2 HQET on the lattice 

In this section we define the effective field theory for QCD containing a heavy quark flavour 
in lattice regularization, starting from the formal 1/m-expansion of the classical theory. 
We drop all terms involving the heavy anti-quark fields as they can be incorporated in 
complete analogy to those containing the heavy quark field iph which we discuss in detail. 1 
Renormalization properties are addressed but the proper choice of renormalization conditions 
and the physics content of the theory is deferred to the next section. 

2.1 Definition of the effective theory 

We consider QCD on the lattice. The explicit form of the gauge field and light fermion action, 
SVei, is not needed for our general discussion, but for some of the following statements to 
hold, an 0(a) improved formulation is required, e.g. the one described in [9]. 2 We denote 
the set of (bare) parameters of the theory with Nf relativistic quarks by Cjv f - Apart from 
the gauge coupling, go, and the quark masses, it will in general also cover some improvement 
coefficients [9]. 

As has been explained by Eichten and Hill [3,4,10], an effective field theory for hadrons 
(at rest) containing Nf — 1 light quarks and one heavy quark (b-quark) with mass m may 
be obtained by a formal 1/m-expansion of the continuum QCD action and the fields, which 
appear in the correlation function under study. The action of the heavy quark is written in 
terms of the four-component field tp^ satisfying 



For simplicity we drop higher-dimensional operators in the effective field theory which involve only light 
quark fields and the gluon field. These terms contribute at higher order in 1/m. 

2 0(o) improvement means that the continuum limit is reached with corrections of 0(a 2 ). 
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Including terms up to the order l/m n in the expansion, the action, discretized on a Euclidean 
lattice, reads 

n 

Shqet = a 4 ^{£ stat (x) + ]T£W(x)}, (2.2) 

£stat(z) = Mx)[V* + 5m]Mx), (2.3) 
&\x) = (2.4) 

where denotes the backward lattice derivative, 5m has mass-dimension one, and the local 
composite fields have mass-dimension 4 + v. Indeed, this form is suggested by a formal 
1/m-expansion at the classical level which yields 3 

5m = 0, (2.5) 
4 1} =? h (-^-B)Vh, ^ 1} = l/m, (2.6) 

4 1} =^hf- ^D 2 )Vh, w^l/m (2.7) 



up to and including the order 1/m. Here, B is a discretized version of the chromomagnetic 
field strength and D 2 a lattice version of the covariant Laplacian in three dimensions. Note 
that a term mip^x^^x) has been removed from the action, since it only corresponds to 
a universal energy shift of all states containing a heavy quark. Removing it makes explicit 
that the dynamics of heavy-light systems is independent of the scale m at lowest order of 
1/m. 

While the action is sufficient to obtain energy levels, for many applications one is inter- 
ested in (e.g. electroweak transition matrix elements) it becomes necessary to also discuss 
correlation functions of composite fields. As an example we take the time component of the 
axial current. In the effective theory it is defined by an expansion similar to eq. (|2.2j) . 



A^ ET (x) = ^A^(x), (2.8) 

A^(x) = a^Af*\x), Ar t (^)=?i(x) 7 o75V'h(^, (2.9) 

A^\x) = J24 v) A^\x), u>0, (2.10) 

i 

where a light quark field, enters and A± is of dimension 3 + v. One may then study for 
instance the correlator (with {^^Tiftjy = tpj'yoT^jo'ipi) 

C^ ET (x ) = a 3 £ (^ QET (-)(^ QET ) f (0)) • (2.H) 



At the classical level the fields are given by 



i 0) = l, A? = ViT^^^h , c4 1} = l/m. (2.12) 



A short derivation may e.g. be found in [11] 
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In general, i.e. at the quantum level, expectation values are denned by a path integral 

(O) = ^ J D[(p]0[(p] e~( Slcl+SllclET \ (2.13) 

Z = J D[Lp]e' {Slcl+SH ^ T) , (2.14) 

over all fields with the standard measure, denoted here by -D [</?]• An important ingredient 
in the definition of the effective field theory is that it is understood throughout that the 
integrand of the path integral is expanded in a power series in 1/m, with power counting 
according to 

up = O (1/m") , af ] = O (1/m") . (2.15) 
In other words one replaces 

exp {-(5 re i + 5'hqet)} = exp I - (s reX + a 4 ^ £ stat (a;)) j x 

{ 1 - a 4 £ (x) + \ [ a 4 £ (*) ] 2 " « 4 E ^ (*) + -} ( 2 - 16 ) 

^ X X X ' 

in eq. ()2.13j) . The 1/m-terms then appear only as insertions of local operators o\ (x) and 

(ar) into correlation functions, and the true path integral average is taken with respect 
to the action in the static approximation for the heavy quark, S = S re \ + a 4 J2 X ^stat(x). 

Power counting leads us to expect that the static theory is renormalizable, requiring 
a finite number of parameters to be fixed to obtain a continuum limit. Indeed explicit 
perturbative [5, 12-15] as well as non-perturbative [16] computations support that this is a 
genuine property of the static effective theory. Would one keep one of the 1/m-terms in the 
exponent, as it is done in NRQCD, renormalizability would be lost and most of what we are 
concluding in this paper would not be true. 

We are still left to discuss the renormalization of expectation values of the type (|2. 13|) 
after inserting the expansion (|2.16|) . This is just the problem of renormalizing correlation 
functions of local composite operators in the static effective theory. Power counting imme- 
diately leads to the conclusion: once all local operators, whose dimensions do not exceed the 
one of the highest-dimensional operator (i.e. v < n) and which have the proper symmetries, 
are included, their coefficients may be chosen such that all expectation values have a con- 
tinuum limit (see e.g. Ref. [17]). Of course, both the operators Cf J \x) in the action and the 
ones in the effective operators such as Jtf\x) have to be included. One may worry that 
due to the sums over all space-time points in eq. (|2.16|) contact terms appear, which lead 
to additional singularities. However, just like in the case of 0(a) improvement discussed 
thoroughly in [9] , the terms needed to remove these singularities are already present once all 
local operators with the appropriate dimensions are included. 

The effective theory is now defined in terms of the set of parameters, 

C m ET = {ck} = C Nt . 1 U{5m}U{4 u) }^W ( f ) }U... , c x = gl. (2.17) 

The ellipses allow for coefficients of further composite operators which will be needed when 
their correlation functions are considered. For the continuum limit of this effective theory to 
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exist, the parameters {ck,k > 1} have to be chosen properly as a function of Note that 
in the notation used here, the renormalization of the effective composite fields is included in 
the set of generalized coupling constants, Chqet- E.g. at lowest order in 1/m, the coefficient 
ckq = Zj^ at is the renormalization constant of the static axial current [14]. 
A few more remarks are in order. 

• Once the proper degrees of freedom, namely the field ip^, have been identified, the 
terms in the effective field theory are organized just by their mass-dimension. The 
expectation that the effective field theory has a continuum limit (is non-perturbatively 
renormalizable) is thus nothing but the usual expectation that composite operators 
mix only with operators of the same and lower dimension. 

• The same argumentation is also the basis of Symanzik's discussion of cutoff effects of 
lattice theories and their removal order by order in a [9,18-20]. An important conse- 
quence is that in general the 1 / m-expansion and the a-expansion are not independent 
but have to be considered as one expansion in terms of the dimension of the local 
operators. If we imagine to start with a theory with a set of operators identified by the 
formal continuum 1 /m-expansion, these operators will for instance mix under renor- 
malization with operators of the same and lower dimension, which are allowed by the 
lattice symmetries but not by the continuum symmetries and which would therefore 
not be in the set of the operators one started with. To avoid this, one has to start 
immediately with the full set of operators of a given dimension, restricted only by the 
lattice symmetries. In other words we have to count a = 0(1 /m). This means also 
that S'rei has to be 0(a) improved to go to order 1/m. 4 

• Of course, symmetries restrict the terms that have to be taken into account. In general, 
out of the space-time symmetries we only have the 3-dimensional cubic group instead 
of the 4-dimensional hyper cubic group. At the lowest order in 1/m there are additional 
symmetries: heavy quark spin-symmetry [21] and the local conservation of heavy quark 
number, which simplify 0(a) improvement (see Sect. 2.2 of [14]). 

• Furthermore it is convenient to formulate the effective theory only on-shell, i.e. for 
low energies as well as for correlation functions at physical separations. Then the 
argumentation of [9] can be taken over literally to show that the equations of motion 
(derived from the lowest-order action) can be used to reduce the set of operators 
C ( i u \x),A ( i u) (x), .... Following the same reference, operators obtained by multiplying 
those of dimension d by a light quark mass are to be counted as separate operators of 
dimension d + 1. 

• Finally note that after using eq. (|2.16j) . the determinant arising from the static quark 
action is just an irrelevant constant. In principle, loop effects of the heavy quark are 
still present in the coefficients c&. 

4 It may be possible to go to higher order in a than in 1/m, when symmetries restrict the allowed mixings. 
An example is provided by O(a) improvement of the static effective theory [14]. Ways to extend this to higher 
orders in 1/m probably exist but we have not investigated this question systematically. 
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2.2 Power divergences and non-perturbative renormalization 

The mixing of operators differing in dimensions by p translates into coefficients diverging 
(when p > 0) as a~ p . In the present context it actually has been checked in perturbation 
theory that these mixings are not forbidden by some accidental symmetry [5]. Due to such 
power divergences, perturbation theory is not sufficient to determine the coefficients c/-. An 
estimate of order would leave a perturbative remainder 

Ac k ~gf +1) a- p ~ a~ p [ln(aA)]-( /+1 ) a -=^ oo (2.18) 

with A the QCD A-parameter. This means that the continuum limit does not exist if the 
coefficients are determined only perturbatively. 

Hence we conclude that a non-perturbative method is needed to determine (at least some 
of) the parameters {cfc}. Such a method will be introduced in the following two sections. 

3 Matching of HQET and QCD 

By QCD we denote the theory including a relativistic heavy quark, the b-quark, while with 
HQET we mean the theory where this quark is incorporated with the action defined in the 
previous section. The latter is an approximation to QCD when the coefficients Chqet = { c fc} 
are chosen correctly. Then we expect 

$ H Q ET (M) = $ QCD (M) + O (1/M n+1 ) (3.1) 

for properly chosen observables, <J>Q CD , in QCD and their counterparts, <£ H Q ET , in the ef- 
fective theory. Amongst the many dependencies of $Q CD we have indicated only the one on 
the heavy quark mass. To be free of any renormalization scheme dependence, we choose the 
renormalization group invariant (RGI) quark mass denoted by M [22]. In order for eq. (|3.1j) 
to hold, all other scales appearing in cj>Q CD are assumed to be small compared to M. Choos- 
ing as a typical low-energy reference scale of QCD the energy scale r^ 1 ^ 400 MeV) [23], 
defined in terms of the QCD force between static quarks, the combination r$M has to be 
large. Thus the symbol 0(1/M n ) is a short hand for 0(l/[roM] n ). 

To give a simple example for a quantity $Q CD , one could take $Q CD = Caa 5 where 

C AA (x ) = Z 2 A a 3 ]T (A)(x)(4))t(0)) (3-2) 

X 

with the heavy-light axial current in QCD, = tpi^^j^iph, and Z\ ensuring the natural 
normalization of the current consistent with current algebra [24,25]. Then eq. (|3.1|) is valid for 
$hqet = e - m z 0C HQET^ with ^hqet^ from eq and in the region l/x < M. 

With the latter, kinematical, condition one takes care that the correlation functions are 
dominated by states with energies (the heavy quark mass being subtracted) small compared 
to M. Furthermore, the factor e~ mx ° accounts for the mass term that had been removed 
from the effective theory Lagrangian as already mentioned after eq. (|2.2[) . Which mass m is 
to be taken here, depends on the convention used to define 5m. As will be explained further 
in Sect. 15.1.11 at each order in l/m, the combination m + ^ln(l + a5m) is uniquely fixed 
by the matching of HQET and QCD. We emphasize again that the same mass m enters all 
correlation functions involving one heavy quark. 



6 



Let us now come to the main problem: determining the parameters in the effective 
theory such that this equivalence between HQET and QCD is true. First of all assume that 
the parameters of QCD have been fixed by requiring a set of observables, e.g. a set of hadron 
masses, to agree with experiment. It is then sufficient to impose 

$ HQET (M) = (J) QCD (M); fc = 1) ... )iVn , (3.3) 

to determine all parameters {ck ,k = 1, . . . , N n } in the effective theory. Observables used 
originally to fix the parameters of QCD may be amongst these <3?^ CD . The matching condi- 
tions, eq. (|3.3|h define the set {ck} for any value of the lattice spacing (precisely speaking, 
for any value of a/ro). 

In principle, each <I> fc could be determined from a physical, experimentally accessible 
observable. However, this would reduce the predictive power of the effective theory since it 
contains more parameters than QCD. Particularly for increasing the order n of the 1/m- 
expansion we then would need to use more and more experimental observables. 

To preserve the predictability of the theory, we may instead insert some quantities 
<3?k CD (M) computed in the continuum limit of lattice QCD. This of course demands to treat 
the heavy quark as a relativistic particle on the lattice, seemingly in contradiction to the 
very reason to consider the effective theory: small enough lattice spacings to do this are very 
difficult to reach. An additional ingredient is thus necessary to make the idea practicable. 
It will be explained in the following section. At this stage the important point is that there 
are no theoretical obstacles to a non-perturbative matching. We end this section with some 
comments on details of the general matching procedure. 

• The observables <J> are assumed to be renormalized. Eq. (|3.3|) is, however, used to fix 
the bare parameters in the action — for each value of a 2 . 

• When one increases the order n in the expansion, new quantities have to be added, 
and at the same time, the parameters of the lower-order Lagrangian, Cj, i < -/V n _i, 
will change in general. This change is due to mixing of the operators and may thus be 
sizeable. 

• It is convenient to take the continuum limit 5 of &® CD before imposing eq. 1)3. 3 j) . If one 
decides not to do this, the lattice spacings on both sides of eq. (|3.3j) should be scaled 
together in order to reach the continuum limit in the effective theory. 

• As mentioned already in the previous section, the terms necessary for Symanzik im- 
provement are taken into account automatically, namely some of the equations (|3,3[) 
may be interpreted as improvement conditions. Working up to the order n, the result- 
ing lattice HQET is correct up to 

error terms = O ((l/m) n+1 ) = O (M^ n+1 \aM) k ^ , k = 0, 1, . . . , n + 1 . (3.4) 

Higher-order terms in 1/M have parametrically larger lattice spacing errors. For ex- 
ample, a treatment of the theory including the next-to-leading operators will give us 
the (l/M)°-terms with 0(a 2 ) errors and the linear 1/M-corrections with 0(a) uncer- 
tainties. Additional work would be necessary to suppress the discretization effects in 
the 1/M-terms to 0(a 2 ). 
5 Or alternatively, work in a sufficiently improved lattice theory and at a small value of the lattice spacing. 
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• There is a close analogy of our proposed matching procedure to what is done when 
the low-energy constants of the chiral effective Lagrangian [26] are determined using 
lattice QCD. An important difference is, however, that the chiral expansion can be 
worked out analytically while here we still have to evaluate the resulting theory by 
Monte Carlo. The reason is that strong interactions remain; the lowest-order theory, 
the static approximation, is non-trivial. 

4 The role of finite volume 

From the theoretical point of view, the matching described in the previous section is sufficient. 
However, we should take into consideration what can be done in a numerical computation. 
To give a concrete example, let us assume that 

(L/af x T/a < 32 3 x 64 (4.1) 

lattices can be simulated, numbers which are realistic for present computations in the 
quenched approximation, but too large for full QCD. We further assume that we deal with 
quantities which have negligible finite size effects when 

L > 2 fm . (4.2) 

Then the smallest lattice spacing reachable is a ~ 0.06 fm, and this number will not be very 
different if the above assumptions are modified within reasonable limits. While such a lattice 
spacing is small enough to perform computations for charm quarks [27,28], the subtracted 
bare mass of the b-quark is about am q ~ 1. In this situation lattice artifacts are expected 
to be very large and it is impossible to obtain the r.h.s. of eq. (j3.3j) . 

The situation becomes quite different when one considers observables defined in finite 
volume with L considerably smaller than 2fm and uses the generally accepted — and also 
much tested — assumption that both QCD and HQET are applicable in a finite volume and 
the parameters in the Lagrangians are independent of the volume. 

4.1 Matching in finite volume 

Instead of eq. (|3.3j) we now consider (remember that N n is the number of parameters in the 
effective theory): 

$HQET (L)M) = $ QCD (L)M)j k =l,...,N n . (4.3) 

This will allow us to have much smaller lattice spacings on the r.h.s. in order to eventually 
approach the continuum limit. A typical choice is L = Lq ~ 0.2 fm. As has been shown in 
the preliminary report of our work [6] , eq. Q4.3|) can be evaluated very precisely for suitably 
selected quantities D and the continuum limit can actually be taken. 

Concerning the l.h.s., we have to take into account that a— effects will certainly be 
significant when the resolution a/I of the finite space-time is too coarse. Hence the lattice 
spacings where the bare parameters {ck(go)} can be determined are a = 0(0.02 fm) = 
O(Lq/10). For such values of a, the computation of the physical observables in the infinite- 
volume theory {L ~ 2fm in practice) would again be impracticable, because lattices with 
too many points (L/a) 4 would be required. Therefore, a further step is necessary to make 
larger lattice spacings and thereby larger physical volumes available in the effective theory. 
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4.2 Finite-size scaling 



Also for this step a well-defined procedure is easily found. First assume that all observables 
3> fe [L, M) have been made dimensionless by multiplication with appropriate powers of 
L. Next we define step scaling functions [8], Fj~, by 

^ QET (sL,M) = F k ({^ QET (L,M), j = l,...,N n \) , k = l,...,N n , (4.4) 



where typically one uses scale changes of s = 2. These dimensionless functions describe the 
change of the complete set of observables {<3? fe } under a scaling of L — ► sL, and we briefly 
sketch how they can be computed. One selects a lattice with a certain resolution a/L. The 
specification of #? T (L, M), j = 1, ... ,N n , then fixes all (bare) parameters of the theory. 
The l.h.s. of eq. 1)4. 4j) is now computed, keeping the bare parameters fixed while changing 
L/a — > L'/a = sL/a. Repeating this for a few values of a/L, the continuum limit of F% can 
be obtained by an extrapolation a/L —* 0. 

HflFT 

An important practical detail is to choose the various quantities such that each 

.Ffc depends only on a few 3> ■ and the bare parameters can be determined rather in- 
dependently from each other. For instance, it is natural to identify the running Schrodinger 
functional coupling g 2 (L) [8,29] with $^ ET (L, M) and to keep all of the light quark masses 
zero in these steps. In this way Qq and the (light) bare quark masses are fixed independently 
of the parameters 5m,io^ u \ af \ . . . coming from the heavy sector. In the quenched approx- 
imation or with two dynamical quarks, the set of bare parameters specifying the relativistic 
sector, CjVf-i, can then be taken over from [22,29-31] without change. 

A few steps — may be two — are necessary to reach a value of L = 0(1 fm), where at the 
same time contact can be made with resolutions a/L that are affordable to accommodate the 
suitable observables on a physically large lattice to realize the original matching condition, 
eq. (|M.3j) . (In our first example, Sect. 15. 1| this role will be played by the B-meson with its 
mass as the physical input.) We note that in principle the size of Lq is rather arbitrary, 
but the following consideration is important. We are matching at a finite value of 1/m and 
a finite order n. Thus the final results will depend on which quantities have been used to 
perform the matching. If one chooses quantities with kinematics where the 1/m-expansion 
is not accurate (or even not applicable) , this will translate into badly determined parameters 
in the effective Lagrangian and large final truncation errors. For this reason, Lq has to be 
chosen such that the 1/m-expansion is applicable which means 

1/L «m, (4.5) 

and Lq cannot be too small. From these considerations it appears that Lq ~ 0.2 fm — 0.4 fm 
is a good choice. 



4.3 Evaluation of the physical observables in the effective theory 

Physical observables usually have to be computed in large volume which, for practical rea- 
sons, means at lattice spacings around 1/20 fm to 1/10 fm. In this region the bare parameters 
of the effective theory are determined as follows. 
One chooses a suitable K such that 

L K = s x L «lfm. (4.6) 
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Iterated applications of eq. Q4.4|) give rise to recursion relations, the solutions of which 
determine quantities Vk = $u (Lk,M) in the larger volume of extent Lx- Next, re- 

l-TOF"T 

garding & k (Lk, M) = 14 as a requirement while setting the number of lattice points to 
Lk/& = O(10) just fixes the bare parameters Chqet- These bare parameters are then known 
at values of the lattice spacing, where the computation of correlation functions in large vol- 
ume is possible in the effective theory and masses and matrix elements can be extracted from 
their large-time behaviour. 

Note that in the notation used here also the renormalization constants of the composite 
operators appearing in the correlation functions are amongst the "bare parameters". All 
quantities are thus renormalized entirely non-perturbatively. 6 One may still wonder how M 
itself is fixed. The answer to this question is provided by the first of the two examples, which 
we will use now to illustrate the general strategy. 



5 Examples 

In this section we supply two applications of our non-perturbative matching strategy of 
HQET and QCD that up to now was formulated in rather general terms: a full calculation of 
the b-quark mass in combined static and quenched approximations fSect. lS"Tj) and a proposal 
for a non-perturbative determination of multiplicatively renormalized matrix elements of the 
static-light axial current, which is different in spirit from Ref. [16] and still awaits a numerical 
investigation. 



5.1 The b-quark mass at lowest order 

Several determinations of the mass of the b-quark, which use the static approximation on 
the lattice (HQET to order (1/m) ), have been published [32]. They all rely on a pertur- 
bative estimate of 5m [33-35] and suffer from a power-law divergence due to the mixing of 
ipfrDoifih and ip^h as discussed in Sect. 12.21 Their precision is thus limited by the fact that 
a continuum limit can not be taken, and it is difficult to estimate the associated uncertainty. 
We here explain how a entirely non-perturbative computation can be done and will also give 
a first result, which can easily be improved in precision in the near future. This step is also 
a prerequisite to perform the matching of other quantities such as the axial current, since 
generically the quark mass enters in the matching step ()4.3j) for all 



5.1.1 Strategy and basic formula 

As indicated already in Sect.|lJ given a resolution a/L, we fix g$ such that the finite- volume 
running coupling of Ref. [29] takes a certain value. Furthermore we set the light quark masses 
to zero (with one exception which will be discussed). In the language of Sect. 0] we have 

^jjHQET = _ 2(L)) (51) 

-HQE 
fc+1 



*ffiT = ^ CAC = 0, fc = l iVi-l, (5.2) 



in terms of the PC AC masses of the light flavour number k, mF CAC , and the running coupling 



g (L) in the Schrodinger functional (SF) scheme [29]. Other choices are possible, but the 

6 This represents an advantage in comparison to [16] where a last step using perturbation theory was 
necessary to get to the "matching scheme" [7, 16], which here we achieve by virtue of eq. 13.31 . 
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above is convenient in view of the present numerical knowledge [22,30,31]. The box length 
L is then parametrized through g 2 (L). A very useful feature of eqs. 1)5. 1J) and (|5.2[) is that 
they do not involve the heavy field at all and determine the bare coupling and quark masses 
independently of the heavy sector; in particular these conditions are independent of the order 
n of the expansion. 

In this section we are only concerned with energies and remain at lowest order in 1/m. 
The only additional parameter in the Lagrangian to be fixed is a5m, i.e. one more condition 
corresponding to k = iVf + 1 in eq. ()4.3j) is needed. We start from a time-slice correlation 
function projected onto spatial momentum zero containing one heavy quark (such as Caa ; 
eq. ()3.2)0 . Denoting it generically as C(xq), in the logarithmic derivative 

r = — In \c(x -a)/C(x + a)] (x /L fixed) (5.3) 
2a L . 

all multiplicative renormalization factors of C{x$) cancel. Below we shall use xq/L = 1/2, 
but other choices are possible. Replacing the fields in the correlation function C{xq) by 
the corresponding effective fields defines C H( ^ et (j;o) i n the effective theory. In the static 
approximation, its logarithmic derivative, r stat , built as in eq. 1)5. 3 J) , depends on 5m in the 
simple form 

-Tstat — -Tstat 

\ Sm =0 + _ m (l + a $ m ) J ( 5 ' 4 ) 

as is easily seen from the explicit form of the static quark propagator. In large volume, which 
due to xq/L = constant also means large Euclidean time, T = T(L,M) will turn into the 
mass of a b-hadron, say m-g. It is now obvious that 

*$*i(L,M) = LT , (5.5) 
$^(L,M) = L(r stat + m) (5.6) 

= L (r stat | 5m=0 + m bare ) , m b are = m + - ln(l + a5m) , 

are sensible assignments to fix the combination mb are via requiring 

M) = *^(L,M) . (5.7) 

Since 5m and m always appear in the combination mbare! they may not be fixed separately, 
unless one arbitrarily defines 5m by an additional condition. 
Due to eq. (|5.4|) . the step scaling function 

a m (f (L)) = 2L [ T stat (2L, M) - r stat (L, M) ) (5.8) 

is independent of mb ar e an d therefore also independent of M; at lowest order in 1/m, energy 
differences in the effective theory do not depend on the heavy quark mass. The step scaling 
function (|5.8j) is thus a particularly simple realization of eq. (|4.4|) . Together with the one for 
the running coupling [8,29], 

a(u) = g 2 (2L)\ s2(L)=u , (5.9) 



it defines the sequence 



= 9 2 (L ) , w = LT stat \ L=Lo , (5.10) 

Uk+l = O-(lifc) > w k+l = + <Tm(Uk) , (5-11) 



11 



which is easily seen to relate F^^Lk, M), with Lk = ^ K Lq, to r sta t(A),M) when the 
sequence no, • • • , uk-i is known: 

K-l 

L T stat {L K , M) = L r stat (L , M) + 2' {k+1) a m (u k ) . (5.12) 

k=0 

Suitable choices for uq and K then allow to arrive at Lk = 0(1 fm). 

Finally one considers the energy E stat of a B-meson in static approximation, given for 
example by 

Ca? ET (^o) xo ~°° AeM-xoE stat ) (L large) . (5.13) 
The energy difference 

AE = E stat - T stat (L K , M) (5.14) 

can be computed with one and the same lattice spacing (i.e. at the same bare parameters) 
for the two different terms on the r.h.s., but of course with different L. Combining eq. (|5,7[) 
imposed in small volume (L = Lq) with eqs. (|5.12j) and (|5.14l) to eliminate Wbare i n m B = 
E s t&t + m bnre (holding in the large- volume limit), we arrive at the basic equation 

A'-l 

L m B = L r(L , M b ) + 2-( fc+1 V m (u fc ) + L AE . (5.15) 

It relates the mass of the B-meson to a quantity T(Lq, M\>), computable in lattice QCD 
with a relativistic b-quark, and the energy differences AE and <7 m which are both defined 
and computable in the effective theory. All quantities on the r.h.s. may be evaluated in 
the continuum limit. Note that all of the terms in eq. ([5.15ft are independent of mbare 
(because the unknown m-bare, and thereby 5m too, drop out in the differences), although 
logically eq. (|5.7|) has been used to fix it non-perturbatively. Our strategy has also been 
presented in a somewhat different way, which is closer in spirit and notation to standard 
HQET applications (but less rigorous), in [36]. 

Eq. (|5. 15f) may be looked at in two different ways. Given the RGI mass of the b-quark, 
Mb, eq. (|5.15|) provides a way to compute the mass of the B-meson. It is more interesting to 
turn this around: taking m-Q from experiment and evaluating (in lattice QCD) T(Lo,M) as 
a function of M, this equation may be solved for -A<fb> Implicitly the bare parameter m-bare 
is thus fixed non-perturbatively, and the problem of a power-law divergence is solved. 

We now give an example for a precise definition of the correlation function C{xq) and 
use the quenched approximation to demonstrate that the continuum limit can be reached in 
all steps while still a very interesting precision is attainable. The reader who is not interested 
in the numerical details may directly continue with Sect. 15.21 

5.1.2 Correlation functions, 0(a) improvement and spin- symmetry 

In our numerical implementation we choose SF boundary conditions with all details as in [16], 
including = 1/2, T = L and C = C = in the notation of that paper. This means that 
0(a) improvement [9, 14] is fully implemented, except for uncertainties in the coefficients 
c^ at , ct and ct originating from their only perturbative estimation. As in [16] it has been 
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checked that the influence of these uncertainties on the observables considered here can be 
neglected compared to our statistical errors. We will therefore not mention 0(a) terms any 
more and perform continuum extrapolations modelling the a-effects as 0(a 2 ). 
For our definition of F we consider the two correlation functions 

fA(xo) = -yE<(^)o(*Kb(y)75Ci(z)}, (5-16) 
kv(xo) = -y E <(H)*(x)? b (y)7fcCi(»)> , (5-17) 

where the label "I" on the axial and vector currents reminds us that their 0(a) improved 
forms are used: 



(Ai)^{x) = ^\{x)^^i) h {x)+acA\{d^ + d%){i) l (x)^'4} h {x)} , (5.18) 
(Vi)fi(x) = tpi{x)j^(x) +acv~(d 1/ + &t) {iip^a^^x)} . (5.19) 



As a consequence of the heavy quark spin-symmetry, their partners in the effective theory 
coincide exactly at lowest order in 1/m and we thus define only one 



a 6 



C hqet (xo) s /r(xo) = __ £ ((Ar%(x) C h (y)75Ci(z)> . (5.20) 



y 5 z 



In these definitions, ^ z Cb(y)7sCi( z ) anci S y z Ch(y)7sCi( z ) are interpolating fields localized 
at the xq = boundary of the SF, which create a state with the quantum numbers of a B- 
meson with momentum p = 0, and for 75 — ► 7/% we have the quantum numbers of a B*. The 
correlation functions are schematically depicted in Fig. [21 more details can be found in [14]. 




Figure 2: Illustration of the correlation function /a- For kv, the insertion of (Ai)o is replaced by (Vi)fe, 
whereas in case of /A at the operator in the bulk is (Af at )o connected to £ h (instead of f b ) by a static quark 
propagator. 

Inserting C{xq) = Ja(xo) arid C{xq) = ky(xo) into eq. ()5.3|) defines Tpg and Fy, respec- 
tively. Their partner in the effective theory is denoted as r s t a t- F)ue to the spin-symmetry, 
either LFps = -L(r sta t + m) or LFy = L(F stat + m) are possible matching conditions at 
lowest order in 1/m. At first order, one has to consider also separately the vector and the 
axial vector correlators in the effective theory since these are split by the a ■ B-term in the 
effective Lagrangian. It is then convenient to define the matching condition as eq. ()5.7|) with 
the spin-average 

r = - A (r PS + 3r v ) . (5.21) 
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With this definition the matching condition (|5.7|) is independent of the coefficient of the 
a ■ B-term also at first order in 1/m, and hopefully 1/m-effects are thereby minimized. 

Having now completed our definition of T and r s t a t , we remind the reader that the mass 
of the light quark is set to zero and thus T is only a function of the heavy quark mass, M, 
and the linear extent of the SF- volume, L (if a-effects are neglected for the moment). 

5.1.3 Matching 

The essential steps of our strategy explained in Sect. 15.1.11 are the matching at L = L$ 
as the starting point, then connecting to L = Lk and from there to the (infinite-volume) 
meson mass. In practice, proceeding from any choice of uq the sequence Uk is only known 
with a certain numerical precision and this has to be taken into account in the error anal- 
ysis. Furthermore we want to take advantage of the numerical results of [22] for triples of 
{L/a, (3, k) corresponding to fixed renormalized coupling and vanishing light quark mass in 
the quenched approximation, as well as of the known function [ro/a](P) and the value of 
imaxAo = 0.718(16) [37], where g 2 (L 

max) — 3.48. So it is convenient to have K — 2 in pre- 
vious formulae and to define (exactly) L2 = 2L max = 1.436 tq. Then, applying the inverse of 
the step scaling function of the coupling [22] twice, we arrive at 

u = 2.455(28) , a~ l {uo) = 1.918(20), (5.22) 

and end up with L = Lq = L2/4 = 0.359 ro ~ 0.18 fm for the linear extent of the matching 
volume. The matching of HQET and QCD is then supposed to be done at u = uq 2.4. We 
could thus keep e.g. u = 2.4484 fixed, where the triples (L/a, (5, k) are known [22], but only 
for 6 < L/a < 16. The spacing of these lattices is still too large to comfortably accommodate 
a propagating b-quark. Instead it is better to work at a constant value of u = 1.8811, varying 
6 < L/(2a) < 16. Nevertheless, T is computed on the lattices with L/a points per direction, 
and the slight mismatch of <t _1 (uo) and 1.8811 will eventually be taken into account together 
with the overall error analysis. 

We now have to determine T(L,M), where in case of the relativistic theory L is always 
to be identified with the extent of the matching volume, Lq, from now on. In order to 
approach its continuum limit, we define 

Q(u,z,a/L) = LT(L,M)\- g2{L/2)=UjLM=z (5.23) 

and extrapolate it as a function of a/L, viz. 

uj{u,z)= lim £l(u,z,a/L), (5-24) 

for a few selected values of z and at fixed u. This requires to compute 0, with z and u 
fixed while changing L/a and therefore also (3. A particular aspect in this step is that in 
imposing the condition of fixed z (at variable /?), the relation between the bare quark mass, 
m q , and the RGI one, M, is needed, where several renormalization factors and improvement 
coefficients enter. Although they had already been determined [22,38], it turned out that it 
was desirable to improve their precision and to determine them directly in the range of (3 
where they are needed in the present application. For this reason they were redetermined in 
Ref. [39], and also w as a function of z, eq. (|5.24j) . was obtained by extrapolation in a/L — > 
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in that work. For the reader's convenience, we reproduce from [39] a graph of the continuum 
values cj(1.8811, z) together with the fit function 

uj(1.8811,z) =a z + ai + a 2 - , a = 0.581 , 01 = 1.226, a 2 = -0.358 (5.25) 

z 

in Fig. |31 In the interval 5.2 < z < 6.6, which is the relevant z-range to extract the RGI b- 
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Figure 3: Continuum limit values of LoT(Lo, M) at fixed coupling g 2 (Lo/2) = 1.8811 as a function of 
z — LqM and its fit function, determined in relativistic QCD but small volume [39]. The b-quark mass scale 
lies near z ss 6. 

quark mass later, this parametrization describes w(1.8811, z) with a precision of about 0.5%. 
A further global uncertainty of 0.9% has to be attributed to the argument z of the function 
co (see Ref. [39]). In order to also take the small statistical error and mismatch in uq into 
consideration at the end, we also need a numerical value for the derivative of u/(1.8811, z) 
w.r.t. u. It was found to be constant in the interesting region [39]: 

= 0.70(1), 6.0<z<6.6. (5.26) 

u=1.8811 

For completeness we also quote the fit result for u;(1.8811, z), if instead of the spin- 
average ()5.21j) it is defined as the continuum limit of the effective energy LTps(L, M). With 
the same fit ansatz as in eq. (|5.25j) . the coefficients then read 

a = 0.587 , ai = 1.121 , a 2 = -1.306 (for T = T PS ) . (5.27) 

The significantly larger a2~term in this case compared to eq. (j5.25|) indicates that the spin- 
averaged combination T has smaller 1/m-errors and should therefore be preferred in the 
implementation of the matching step. 7 

5.1.4 Finite-size scaling 

The next step is the numerical determination of the step scaling function <r m , eq. (|5.8|1 . 
and then of the (M-independent difference) Lq [r stat (Lk , M) — r stat (Lo, M)] as given by 
eq. (|5~T2l . 
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Figure 4: Lattice step scaling function E m and its continuum limit extrapolation linear in (a/L) 2 , which 
uses only the four smallest values of a/L. 



u 


cr m (n) 


2.4484 


-0.205(18) 


2.7700 


-0.133(26) 


3.4800 


0.040(25) 



Table 1: Results for the continuum step scaling function <7 m (it). 



At finite lattice spacing, the step scaling function is defined by 

E m («,a/L) = 2L[r stat (2L,M) -r stat (L,M)]| g2(L)=M , (5.28) 

where, as mentioned earlier, the (light) quark mass is set to zero. The exact definition of the 
massless point does not play an important role; it is as in Ref. [16]. In fact, we evaluated 
S m directly from the correlation functions computed there, where details of the simulations 
can be found, too. Numerical results for E m are listed in Table IA.1I in the appendix. They 
are extrapolated to the continuum limit via 

a 2 

S m (u, a/L) = a m (u) . (5.29) 

Given our data at various resolutions, this is a safe extrapolation (cf. Fig. |1J) leading to 
continuum values reported in Tabled Setting now u as in eq. (|5.22[) . we want to compute 
(recall that K = 2 now) 

L [r stat (L 2 ,M) -r stat (L ,M)] = -a m (u ) + -a m (ui) , (5.30) 

with u\ = cr(uo) = 3.48(5) [22]. To this end it is convenient to represent the data of Tabled 
by a smooth fit function, 

<7m(w) = S + SlU + S2U 2 , (5.31) 
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Figure 5: Continuum step scaling function a m (u) and its polynomial fit. 

shown in Fig. |SJ This fit implies the numerical value 

2<7m(«o) + 4^m(«l) = -0.092(11) (5.32) 

for the combination (|5.30|) . Here the uncertainties in uq,u\ are neglected, since they would 
contribute only a negligible amount to the total error of eq. (j5,32|) . Dropping the S2tt 2 -term 
in eq. (|5.31[) would give indistinguishable results. 

5.1.5 The subtracted B-meson mass 

As a last ingredient for the basic formula in eq. (|5.15|) we have to address LqAE, with AE 
the energy difference between the B-meson's static binding energy and the effective energy 
of the static-light correlator fj£ , eq. (|5,14|) . at the same values of the lattice spacing. We 
evaluate this quantity starting from results for aE stat (go) in the literature [40]. They are 
interpolated in the mass of the light quark to the strange quark mass (see also Sect. 5.2. of 
Ref. [16]) and then correspond to a B s -state. Since O(a) improvement was not employed 
in the computation of Ref. [40], we also need r stat (-L2) for the unimproved theory. Given 
the simulation results reported in Appendix C.2. of Ref. [16], aT stat (gQ, L/a) with L = L2 = 
1.436 ro is straightforwardly obtained for 5.68 < = 6/ ' g$ < 6.5. These numbers are collected 
in Table lA~2l of the appendix and well described by 

ar stat ( 5o , L/a)\ L=L2 = 0.394 - 0.055 (J3 - 6) - 0.218 {(3 - 6) 2 + 0.229 (/3 - 6) 3 (5.33) 

with an absolute uncertainty of less than ±0.002 in the range of (3 mentioned. 

The combination LqAE is shown in Fig. El Its errors are dominated by those of aE stat . 
Since they are rather large and also the lattice spacings are not very small, we refrain from 
forming a continuum extrapolation. Instead we take 

LqAE = 0.46(5) (5.34) 

7 In the preliminary computation of the b-quark's mass reported in refs. [6,7], the matching was performed 
using only the energy in the pseudoscalar channel, Tps. 
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Figure 6: Subtracted, dimensionless B a -meson energy evaluated from the bare numbers of [40]. 

as our present estimate. As seen in the figure, its error encompasses the full range of results 
at finite a and the true continuum number is expected to be covered by it. 

No doubt, a continuum limit with small error (at least a factor 3 smaller) can be achieved 
here in the near future, incorporating 0(a) improvement and using the alternative discretiza- 
tion of static quarks of Ref. [41]. 



5.1.6 Determination of the quark mass 

Now we are in the position to put all pieces together and solve the basic equation (|5.15|) for 
the b-quark mass. This amounts to determine the interception point of the function 



uj (o-- 1 (uq),z) =w(1.8811,z) +A 



8 ( ^ 



(5.35) 



m=1.8811 



where Au = (a 1 (uq) — 1.8811) = 0.037(20) accounts for the slight mismatch in the couplings 
fixed, with the combination 



Wstat = L m B - { \ cr m (u ) + \ <7 m (iti) } - L AE . 



(5.36) 



All of the necessary quantities are known from the foregoing three subsections and the 
experimental spin-averaged B-meson mass me = m B s = 5.40 GeV is taken as the physical 
input. As illustrated in Fig.[7J this procedure yields a value for LqM^ with an error. Together 
with Lo/ro = L mgx /(2ro) = 0.359 it is converted to our central result 



r M h = 16.12(25)(15). 



(5.37) 



Here the second error in parentheses stems from the additional 0.9% uncertainty of z in 
u(u, z) that was mentioned in the context of eq. (|5.25|) . With rg = 0.5 fm this translates into 



M b = 6.36(10)(6) GeV , m b (m b ) = 4.12(7)(4) GeV , 



(5.38) 



where the running mass, m^), is evaluated with the four-loop renormalization group func- 
tions and in the MS scheme. 
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Figure 7: Graphical solution of the basic formula, eq. 15.15L for the dimensionless RGI b-quark mass, 
Zb = LoMb- The contributing pieces are repeated explicitly in eqs. ^5.350 and l|5.M6[l . 

One should remember that this result is valid up to 1/m-corrections and it has been 
obtained in the quenched approximation. One may estimate the ambiguity due to setting the 
energy scale in the quenched approximation as usual. Varying the value of ro in fm by 10%, 
we obtain changes of about 260 MeV in and 150 MeV in We emphasize, however, 

that the scale ambiguity can serve only as a rough guide to the impact of quenching. Finally, 
we also want to stress once more that in Sect. 15.1.51 the contribution from the subtracted 
B-meson mass to our result on the b-quark mass is still based on ordinary Wilson fermion 
data and will be hopefully much improved soon. 

5.2 The B-meson decay constant at lowest order 

To lowest order in 1/m, the bare matrix element 



evaluated with the static action together with a renormalization factor Z^ at (gQ, oMb) allows 
to determine the B-meson decay constant in static approximation [3]: 



The renormalization constant Zj£ (go, aM^) has already been computed in the quenched ap- 
proximation [16], using a quite different method to the one described in this paper. However, 
the method of Ref. [16] is not easily extended to include 1/m-corrections. Therefore, as a 
second example of the use of our general strategy, we here describe an alternative method 
which may be extended to include 1/m-corrections. 

The key formula, valid up to corrections of 0(l/m), has already appeared and was 
briefly discussed in Ref. [7]: 




(5.39) 




(5.40) 




X 



$Q CD (L ,M b ). 



(5.41) 



19 



In the rest of this section we give precise definitions of the various factors and explain the 
formula in the general framework of Sects. El and 0] We assume that the b-quark mass is 
already known. Terms necessary for 0(a) improvement are not written explicitly; they can 
easily be added. 

5.2.1 Matching 

As a preparation for the matching of the axial current and the associated determination 
of Z A tat = ckq , we start from f^ at (xo), defined in Sect. 15.1.21 This correlation func- 
tion renormalizes multiplicatively: (/ A tat )R = Z^ h Z^ 1 Z| tat /| tat . In order to eliminate the 
renormalization factors Zr h ,Z& of the boundary quark fields, we consider in addition the 
boundary-to-boundary correlation, 

/f at = -|^ E (Cl'^TsCh'^CU^TsClW) , (5.42) 

u,v,y,z 

which is renormalized as (/i tat )R = (Z^Z^) 2 ff tat . In the ratio 

X ^LIa) S S^m (5.43) 

the unwanted Z-factors cancel, and with the choice xq = L/2 it is also independent of the 
linearly divergent mass counterterm 5m. The renormalized ratio is denoted by 

$HQET (LMb) _ x R (u,z h ,a/L) (5.44) 
= Zf at (g ,aM h ) X(g , L/a)\ s2(L)=u , z h = LM h . 

Here, as wherever we do not mention the light quark masses explicitly, they are set to zero. 

In QCD we define the corresponding quantities (Za denotes the standard renormaliza- 
tion constant for the QCD axial current [24,25] and /i is the analogue of /f tat with £h — > Cb) 

Y(g ,z h ,L/a) = t^fil , (5.45) 
v h 

Y R (u,z h ,a/L) = Z A (g )Y(g ,z h ,L/a)\ § 2( L)=u , (5.46) 

$ QCD (L,M b ) = lim Y R (u,z hl a/L), (5.47) 

a/L^O 

and the matching equation to be imposed in the small volume (of extent Lq) is 

$ HQET (L ,M b ) = $ QCD (L ,M b ) with g 2 (L ) = u = fixed. (5.48) 

While in agreement with the notation in previous sections $ H Q ET on the l.h.s. has a lattice 
spacing dependence that is only implicit (cf. eq. ()5.44j0 . for the r.h.s. the continuum limit is 
to be taken (cf. eq. (|5.47|) ). In the quenched approximation, the particular value for no may 
be chosen as in Sect. 15.11 
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5.2.2 Finite-size scaling 

Computing the step scaling functions built as 

X(g ,2L/a) 



crx(u) = lim 



(5.49) 

9 2 {L)=u 



a/L->o X(g ,L/a) 
allows to reach larger values of L via the recursion 

$HQET (2L)Mb) | a=o = ax ( ^2 (L)) x <|>HQET (L)Mb) | a=o ^ 

We note in passing that in eq. (|5.49f) we could have written X — > Xr, since the same Z| tat 
enters in numerator and denominator. For the following we assume eq. (j5.5U|) to be iterated 
K times, connecting $ HQET (L , M b ) to <& UQET (L K , M b ) with L K = 2 K L , i.e. in a numerical 
implementation, details will be very similar to those described in Sect. 15.1.41 

5.2.3 The decay constant 

To finally arrive at F-Qy/m-Q, one defines the renormalization constant in view of eq. I|5.44[) . 

[90, aM h ) - — — — — — — , u K - g {L K ) , (5.51) 

X{go,L K /a) X(g ,L K /a) 

and — bearing in mind that $ H( ^ ET (Lx, Mb)| a =o via ec i- (|5.5U|) and the matching condition 
(|5.48|) may be evolved back to the renormalized quantity $^ CD (Lo, Mb) in QCD — replaces 
Zf £Lt (g ,aM h ) in eq. (|5.4UI) by rewriting: 

FbV^ = p(uk) x $ HQET (L x ,M b )| a=0 , (5.52) 

M(g Q ) 



p(u) = lim R(u,a/L), R(u,a/L) 



a/L^o ' ' X(g ,L/a) 



9 2 {L)=u 



(5.53) 



Although we have chosen only u and a/L as arguments for R, it does depend on the masses 
of the light quarks used for the evaluation of the matrix element Ai. These masses have to be 
set (or extrapolated) to the physical ones to obtain the correct matrix element in question; 
conveniently, the light quark masses are put to zero in all other quantities as mentioned 
above. In this way the effective theory is renormalized in a (light quark) mass-independent 
renormalization scheme. 

Choosing as an example K = 2, we may combine all ingredients into the expression 

^Bv^b = p(u 2 ) x ox(«i) x ^x(no) x $ QCD (L , M h ) , (5.54) 

where the various factors correspond one to one to those in eq. 1)5. 41|) but are now rigorously 
defined and can be computed in the continuum limit. Our notation is most appropriate for 
the lowest order in 1/m. At higher order, additional matching equations and step scaling 
functions have to be defined and ax will depend on Mb, which is not the case at lowest 
order in 1/m. In fact, at this order all the (heavy quark) mass dependence of F-Qy/m-Q is 
contained in <I>Q CD (Lo, Mb) that is calculable in small- volume lattice QCD with a relativistic 
b-quark [42]. 
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5.2-4 Relation to other approaches 

We finally want to compare the strategy of renormalizing Ag tat , presented in this paper, to 
the one of [14,16], which followed the general ALPHA Collaboration strategy of obtaining the 
renormalization group invariant matrix elements of composite operators non-perturbatively 
and then using (continuum) perturbative results to find the operators normalized in the 
matching scheme at finite renormalization scale. (A recent review is Ref. [7].) For the 
application to the HQET, the natural renormalization scale is then the mass of the b-quark. 

Remaining strictly at lowest order in 1/m, there is no mixing with lower-dimensional 
operators. Consequently, perturbation theory can be applied. In particular the large-mass 
behaviour of i ? By / ^B, which is given by the mass dependence of Zf^ at , is computable leading 
to the asymptotics [43, 44] 8 

$gg< = constant, (5.55) 

1 11 - §(iV f - 1) 
70 b °~ 16^ • 

Furthermore, the function Cps(Mb/A^jg) = ^b^/^b/^rgi ^ s known perturbatively up to 
and including g 4 (mt,)-corrections to the leading-order equation Q5.55|) [45-50]. These per- 
turbative computations, like our non-perturbative method, are based on the renormalization 
of the static axial current where the finite part is determined by matching, eq. (|3.3|) , We 
denote this renormalization scheme by the "matching scheme" [7, 16] . 9 

A finite, renormalized static axial current can of course also be defined by other renor- 
malization conditions involving a renormalization scale /i in a suitably chosen intermediate 
scheme. Matrix elements $i n t er (^) of the renormalized current in this scheme will then not 
necessarily satisfy eq. (|3.1|l . but it is easy to see that 

lrm [8vr6 a( / u) ] _7o/(2f,o) = *asi , « = t 7(4*0 , (5.56) 

is independent of the intermediate renormalization scheme. In Ref. [16] a finite-volume 
scheme was adopted, which allows to evaluate the limit in eq. (|5.56l) through some finite-size 
scaling steps followed by perturbative evolution at very high \x. The results obtained in 
this way are then combined with the perturbative approximation of Cps(Mb/Aj^g). Owing 
to this last step, they are accurate up to relative errors of order cf'im^). This discussion 
should have made evident that the essential difference of the method presented here is not 
the absence of these perturbative errors, which are expected to be quite small, but rather 
the tempting possibility to include 1/m- corrections. 

6 Uncertainties and perspectives 

Following the general strategy introduced in this paper opens the possibility to perform clean 
non-perturbative computations using the lattice regularized HQET. The dangerous power- 
law divergences are subtracted non-perturbatively through the matching in small volume. 

8 The leading-order coefficient bo of the /3-function is taken for Nf — 1 flavours, since the b-quark does not 
contribute. Note that here N{ — 1 = corresponds to the quenched approximation. 

9 Non-perturbatively, the matching scheme is unique up to higher orders in 1/m; in perturbation theory, a 
residual scheme dependence due to the choice of the other renormalized parameters, such as g, remains. 



km [ln(M h /A m )) Wb0) F B ^ 
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This is not only a theoretical proposal: in Sect. 15.11 we showed that in a concrete physics 
application the statistical errors of Monte Carlo results are quite moderate. In fact they can 
be expected to be even smaller, when an alternative discretization of the static approximation 
is employed [41]. 

We emphasize that the result for the b-quark mass in Sect. 15.11 is valid up to 1/Mb- 
corrections. If we had used V = Tps instead of the spin-averaged energy in the matching 
step, a 0.4 GeV higher number for Mb would have been obtained. We do not regard this 
shift as a realistic estimate of the magnitude of 1/Mb -corrections but believe that they 
are significantly smaller, as indicated by the smallness of the associated coefficient 02 in 
the numerically determined quark mass dependence of the spin-averaged energy, eq. (j5.25|) . 
Nevertheless it is clear that a precision determination of Mb requires to take the 1/Mb- 
corrections into account, even if only to show that they are small. 

In general one may argue that the matching step (carried out at order 1/A1£) contains 
l/(Lo-^b) n+1 - uncer t arn ti es m addition to the unavoidable l/(roMb) n+1 -terms (we remind 
the reader that we take 1/ro ~ 0.4 GeV as a typical QCD scale). Whether or not the former 
terms are larger than the latter can only be decided if the linear extent of the matching 
volume, Lq, is varied. Increasing it would demand even smaller lattice spacings. 

Since 1/Mb -corrections are computable, they should be determined to arrive at precision 
computations for B-physics observables, e.g. for the phenomenologically interesting B-B 
mixing amplitude. On the other hand we expect (l/Mb) 2 -corrections to be very difficult 
in practice, because they require also 0(a 2 ) improvement of the whole theory. Fortunately, 
(l/Mb) 2 -corrections do not appear to be very important [42] but, as in every expansion, this 
issue has to be studied case by case. 

An attractive property of our strategy is that it does not involve any particularly large 
lattices and therefore all the steps outlined in the present work can also be performed with 
dynamical fermions. These calculations are presumably no more difficult than for instance 
those of Refs. [51-53] in the light meson sector. 
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Appendix A Detailed numerical results 

In this appendix we collect some of the numerical results underlying the computation of the 
b-quark's mass in Sect. 15.11 

The numerical data on the static-light meson correlator in the Schrodinger functional 
(SF), which are necessary to evaluate its logarithmic derivative r stat , see eqs. (|5.3(l and (|5.2U|) . 
have already been obtained in the context of the non-perturbative renormalization of the 
static-light axial current, Ref. [16]. Hence we refer the reader to this work for more details. 
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-11 T\ 


o nil 

P = 6/So 


K 


L/a 


ar stat (L) 


or stat (2L) 


L m (u,a/L) 


2.4484 


6.7807 


0.134994 


6 


0.3082(6) 


0.3123(6) 


0.049(10) 




7.0197 


0.134639 


8 


0.2899(4) 


0.2881(7) 


-0.028(13) 
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0.134380 


10 
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0.2621(5) 
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16 


0.2444(2) 


0.2396(4) 


—0.151(16) 


2.7700 


6.5512 


0.135327 


6 


0.3282(6) 


0.3374(10) 


0.111(14) 




6.7860 


0.135056 


8 


0.3068(8) 


0.3098(9) 


0.048(19) 




6.9720 


0.134770 


10 


0.2916(4) 


0.2907(6) 


-0.019(15) 




7.1190 


0.134513 


12 


0.2795(4) 


0.2769(5) 


-0.060(16) 




7.3686 


0.134114 


16 


0.2592(6) 


0.2570(8) 


—0.068(32) 


3.4800 


6.2204 


0.135470 


6 


0.3629(6) 


0.3952(7) 


0.388(12) 




6.4527 


0.135543 
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0.3399(5) 


0.3554(11) 


0.249(19) 




6.6350 


0.135340 


10 


0.3219(6) 


0.3311(10) 


0.185(23) 




6.7750 


0.135121 


12 


0.3081(3) 


0.3141(7) 


0.143(20) 




7.0203 


0.134707 


16 


0.2858(3) 


0.2883(7) 


0.080(24) 



Table A.l: Results for the lattice step scaling function E m . For the intermediate quantity r sta t(i) the 
argument M is suppressed, since it is evaluated at rrtbarc = and so does not depend on M. 



In Table HOI we list the numerical results on the lattice step scaling function E m , defined 
through eq. (|5.28p . 

Another ingredient in the determination of is the subtracted B-meson energy AE, 
eq. (|5.14|) . which amounts to calculate the static effective energy r stat at L = 1.436 tq. Our 
quenched results for this quantity, both for the 0(a) improved case (with non-perturbative 
c sw from [54] and the 1-loop value [55] for the coefficient c^ at to improve the static-light 
axial current) as well as for the unimproved theory (where both improvement coefficients are 
set to zero), are shown in Table lA~2l and Fig. IA.ll These numbers are well represented by 
the polynomial expressions 



aT s ta,t(go, L/a)\ L=1A36 



ro 



0.410-0.132 (/3-6) 

for 6.0 < j3 = 6/#q < 6.5 and c sw 

0.394 - 0.055 (13-6)- 0.218 [f3 - 
for 5.68 < (3 = 6/g% < 6.5 and c s 



= non-perturbative 

6) 2 + 0.229 ((3 - 6) 3 
. = 



where their absolute uncertainty is below ±0.001 and ±0.002 in the indicated ranges of 0, 
respectively. In Sect. 15.1.51 we restrict the analysis only to the case of unimproved Wilson 
fermions, c sw = 0, but the 0(a) improved parametrization may be required when also 0(a) 
improved data on the static binding energy -E s tat at various lattice spacings will become 
available. 
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